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Abstract 

After giving the most general formulation to date of the notion of integrability for axially 
symmetric harmonic maps from R 3 into symmetric spaces, we give a complete and rigorous 
proof that, subject to some mild restrictions on the target, all such maps are integrable. Fur- 
thermore, we prove that a variant of the inverse scattering method, called vesture (dressing) 
can always be used to generate new solutions for the harmonic map equations starting from 
any given solution. In particular we show that the problem of finding 7V-solitonic harmonic 
maps into a noncompact Grassmann manifold SU(p,q)/S(U(p) x U(q)) is completely reducible 
via the vesture (dressing) method to a problem in linear algebra which we prove is solvable 
in general. We illustrate this method by explicitly computing a 1-solitonic harmonic map for 
the two cases (p = l,q = 1) and (p = 2, q = 1); and we show that the family of solutions 
obtained in each case contains respectively the Kerr family of solutions to the Einstein vacuum 
equations, and the Kerr-Newman family of solutions to the Einstein-Maxwell equations. 

1 Introduction 

Unlike Maxwell's linear electromagnetic field theory, which is generally solvable by linear superpo- 
sition of plane-wave solutions, the explicit solvability of geometric field theories such as principal 
chiral field models, non-linear sigma models, Yang-Mills connections and Einstein equations, is se- 
riously hampered by their highly nonlinear character, unless the number of independent variables 
can be reduced to 2 via suitable symmetry assumptions. In that case all of the above mentioned 
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geometric field theories are completely integrable [8, 9, 5, 26, 35, 33, 20]. Recall that a nonlinear 
system of PDEs is said to be completely integrable if it can be cast as the compatibility condition(s) 
of an overdetermincd linear system of differential equations, called a "Lax Pair" or Lax System (see 
[31, 3] and references therein). 

Interestingly, all of the geometric field theories named above, in their complete integrability 
regime, are instances of harmonic maps. Such maps generalize the notion of geodesies to higher- 
dimensional domains, and are the simplest of all nonlinear geometric field theories. The Principal 
Chiral Field model involves harmonic maps into Lie groups. Its integrability has been addressed in 
[56, 54, 46]. Nonlinear Sigma-Models are nothing but harmonic maps into symmetric spaces. Their 
integrability has been studied in [38, 56, 18, 40, 50]. More generally, harmonic maps from surfaces 
into compact symmetric spaces and their connections to integrability have also been explored [45, 
44] . On the other hand, stationary axisymmctric Einstein Vacuum and Einstein-Maxwell Equations 
reduce to harmonic maps into the real and complex hyperbolic plane, respectively [33, 48, 49], and 
thus the study of their integrability [8, 9, 5, 20] can also be subsumed into that of harmonic maps. 
Here we recall that by "stationarity" and "axisymmetry" of a solution to the Einstein equations 
of general relativity and gravitation one means the existence of two Killing fields, one timelike 
generating an R action and the other spacelike generating a circle action, for the metric tensor of 
a four-dimensional Lorentzian manifold, which is the unknown in these equations 1 . 

In this paper we show that the integrability results for the above-named models are special 
instances of a more general theorem, which we prove, namely that axisymmetric harmonic maps of 
R 3 into symmetric spaces are completely integrable. 

A key feature of completely integrable systems is the possibility of implementing the Inverse 
Scattering Mechanism (ISM) to find new exact solutions from old ones. This has been exploited in 
particular for scalar equations, such as the Korteweg-de Vries (KdV) hierarchy [24, 2], and the cubic 
nonlinear Schrodinger equation [52, 54], for which the inverse scattering transform is computable by 
solving the Gelfand-Levitan-Marchenko integral equations. These integral equations are however 
less susceptible to explicit treatment when dealing with second-order evolution equations such as 
sine-Gordon [1, 53, 55, 41], or non-scalar- valued systems such as harmonic maps. In the case of sine- 
Gordon, alternatives to ISM such as the Backlund transform, have proved their utility as solution- 
generating mechanisms. Here we will develop a workable approach for obtaining new harmonic maps 
from old ones by supplying the rigorous foundation for a method akin to the Backlund transform, 

1 If both of the two Killing fields are spacelike, then the equations reduce to a wave map, i.e. the hyperbolic 
analogue of a harmonic map. One could also assume other types of actions for the Killing fields, e.g. both of them 
generating circle actions, etc. In the context of globally hyperbolic and asymptotically fiat spacctimes, the most 
natural choice is the one we are considering here. Our results however generalize without difficulty to wave maps 
and other cosmological situations. 
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called the vesture or dressing technique, explained below. 

Expressed informally, we summarize the main results of this paper as follows: 

THEOREM. Let G be a real semisimple Lie group and let K be a maximal compact subgroup 
of G. Then any axially symmetric harmonic map from K 3 into the Riemannian symmetric space 
G/K satisfies an integrable system of equations. Furthermore, it is always possible to generate new 
harmonic maps from a given one using the vesture (dressing) method. 

A precise statement appears in Theorem 5. 

The results of our paper go beyond the current literature by providing a general framework 
for the study of harmonic maps into the noncompact symmetric spaces commonly appearing in 
mathematical physics. We also provide mathematically rigorous proofs, both of the integrability 
claim and also of the solvability of the algebraic equations obtained through the dressing technique, 
under quite general assumptions. To our knowledge, many of the issues we have encountered and 
overcome in providing these rigorous proofs have not been previously addressed in the literature. 

We furthermore demonstrate how the dressing technique can be employed to construct by purely 
algebraic means, a 2nA-parameter family of harmonic maps, for any integer N, into any noncompact 
Grassmann manifold SU(p,q)/S(U(p) x U(q)) with p + q = n, starting from any given harmonic 
map. As an explicit example, this task is then carried out for N = 1 in the two cases (p = 1, q = 1) 
and (p = 2,q = 1). In both of these cases the initial solution, the so-called "seed", when viewed as 
the metric of a stationary axisymmetric spacctime, corresponds to the Minkowski metric. We show 
that in the first case, the four-parameter family of the so-called "1-solitonic" maps one obtains, 
when viewed as a solution of Einstein's Vacuum Equations, contains the Kerr family of spacetime 
metrics [29], while in the second case the six-parameter family obtained contains the Kerr-Newman 
family [36]. We thus provide a complete, rigorous, and at the same time concise derivation of two 
of the most significant exact solutions of Einstein's Field Equations. These examples also suggest 
that the approach via dressing may make the task of generating meaningful solutions for other 
effectively two-dimensional geometric field theories much more tractable. 

We now briefly summarize the main points of inverse scattering and vesture methods in the above 
context. The classical ISM can be described heuristically as follows: given a nonlinear first-order 
evolution equation 

ut = F(u,u x ,u xx , . . .), (1) 

for a scalar function u — u(t,x), one considers an associated eigenvalue problem Lip = \ip for an 
isospectral family of linear differential operators, i.e. L — L(t) such that A t = 0. In the KdV case, for 
example, L = —-§^i + u(t,x). The direct scattering problem consists of finding a scattering matrix 
S(t,X) with the property that (loosely speaking) lim ip(x,t,\) — S(t,X) ■ lim ip(x,t,\). Note 
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that S is a matrix because the asymptotic eigenspaces are multi-dimensional. Given Cauchy data 
u(0, x) for (1) one may use direct scattering to find S(0, A). Now, it turns out that the evolution in 
t of S(t, A) is governed by a linear equation, and also that the eigenfunctions ip satisfy a secondary 
equation -^ip — Bip, where B = B(u,u x ,u xx , . . .) is another differential operator. The isospectral 
condition implies that B satisfies a compatibility condition with L, namely L t = [B,L], which 
agrees precisely with the nonlinear evolution equation (1) of interest. For such cases, by solving 
the inverse scattering problem, i.e. the Gelfand-Levitan-Marchenko (GLM) integral equations, one 
recovers the potential u(t, x) from a given S(t, A), thus solving the original evolution equation (see 
Figure 1). PDEs to which this procedure applies, e.g. KdV, are referred to as being integrable by 
way of ISM. 

u(x,0) > 5(0,A) 
l scattering 

4- 

, inverse , N 

u(x, t) < — S(t, A) 

scattering 

Figure 1: Classical ISM 

This technique may also be utilized to address certain second-order evolution equations, such 
as the sine-Gordon equation. In that case the PDE appears as the compatibility condition 



for a system of matrix equations 



U t -V x + [U, V] = 



for a matrix function ^ = ^(x, t, A), where U, V are given 2x2 matrices depending on the solution 
u(t, x) of the sine-Gordon equation, as well as on the so-called spectral parameter A. The Backlund 
transform, for example, can now be used in place of solving the GLM integral equations, to obtain 
new solutions for the sine-Gordon equation. 

A solution-generating mechanism closely related to the above procedure for sine-Gordon was 
introduced in [8] to treat Einstein's vacuum equations: Assuming existence of two commuting 
Killing fields generating a continuous group G of isometries for a Lorentzian manifold (A4,g), the 
vacuum equations 

= 0, a*, i/ = 0, ... ,3 

where R M „ is the Ricci curvature tensor of g, can be viewed as the compatibility conditions for a 
linear evolution problem 

£* = A*, *| A=0 = g' (2) 
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where g' = g'(x) is the metric of the 2-dimcnsional quotient manifold A4/G, £ is a matrix operator 
involving differentiation in the complex parameter A as well as in x, and A is a matrix depending 
on A and on g'. In [8], the above system (2) was shown to be the Lax system for the reduced 
Einstein's Vacuum equations, which established their complete integrability. 

As mentioned before, complete integrability alone does not suffice to show that there is a work- 
able solution-generating method, because of inherent difficulties in solving the GLM equations. In 
[8] the authors provide an alternative: One first chooses an "initial seed" metric go and solves (by 
any means possible) the linear system (2) to obtain a generating matrix ^>q. The evolution prob- 
lem for the scattering matrix S in the classical approach is now replaced by a dressing or vesture 
technique, in which new solutions g' are constructed by first solving a linear system of algebraic 
equations for the dressing matrix \i which has the property that = x^o is a solution of (2), and 
then setting the parameter A equal to zero in order to recover g' (see Figure 2); the metric g on M. 
solving the original vacuum equations is then recoverable from g' by quadratures. This procedure 
for generating g' will be generalized and fully explained in Section 3. 

£*=A* , T . / / \\ 



*U=o=So 



l 



*=X*o 



g< — *(g',A) 

A=o, quad. 

Figure 2: Vesture for Einstein's Equations 

Incidentally, it was known already [21, 22] (even though not cited in [8]) that the Einstein 
vacuum equations in the stationary, axisymmetric case reduce to a single equation in terms of a 
complex- valued scalar function, called the Ernst potential. It turns out that this is the equation for 
an axisymmetric harmonic map from M 3 into the hyperbolic plane, which is the simplest example of a 
non-compact Grassmann manifold: H R = SL(2,R)/SO(2) = SU(l, 1)/S(U(1) x U(l)). Our thesis 
is that the integrability of reduced Einstein equations is simply a special case of a more general 
phenomenon, namely the integrability of axisymmetric harmonic maps from K 3 into symmetric 
spaces, which is the focus of our paper. 

The rest of this paper is organized as follows: In Section 2 we cover the preliminary background 
needed for the study of harmonic maps into Riemannian symmetric spaces. Section 3 is devoted 
to establishing integrability of such maps, and showing how the vesture method is implemented for 
them. In particular we provide the first rigorous proof that the resulting linear algebraic system is 
solvable in general. In Section 4 we specialize to the case of noncompact Grassmann manifolds G P , q 
and show how the problem of finding N-solitonic maps into them is reduced to solving a 2 N x 2N 
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linear system. We then carry out the computation explicitly for the case when the target of the 
map is either Q\^\ or ^2,1 and N = 1. In each case we choose a starting solution that corresponds 
to the Minkowski metric, and show how to obtain the Kerr, respectively Kerr-Newman solution in 
this way. Avenues of further exploration are briefly discussed at the end of the paper. 

2 Harmonic Maps into Lie Groups and Symmetric Spaces 

In this section, we establish our notation and introduce the necessary terminology for defining 
harmonic maps into Lie groups and symmetric spaces. 

2.1 Lagrangian field theory 

We adopt the approach of [14] for the general set-up: Let (Ai m , g) and (Af n , h) be two Riemannian 
or pseudo-Riemannian manifolds. Any differentiable mapping / : M. M can be viewed as a 
section of the velocity bundle V= [J L(T x M,T q Af), where L(V, W) denotes the space of 

linear transformations from vector space V to vector space W. Using local coordinates (x 7 q,v^) 
on V, such a section is given by sj(x) = (x, f(x) 1 Df{x)). A Lagrangian L is an m-form defined on 
V, i.e. L = £(x, q, v®)e[g] where e[g] is the volume form of (M, g). Once evaluated on a section Sf, 
the Lagrangian becomes an m-form on M, so that it can be integrated on a domain V in M, and 
the resulting functional is called the action corresponding to L: 

A[f,V] = f Losf. (3) 
Jv 

A critical point of the action A, with respect to variations that are compactly supported in V, 
is a solution to the Euler-Lagrange equations for L in V. By analogy with classical Hamiltonian 
mechanics, the quantities (v^) are called the canonical velocities, and their duals with respect to 
the Lagrangian density I are called the canonical momenta 

Pa ' 9< 

The canonical stress is by definition the Legendre transform (with respect to the velocities) of the 
Lagrangian density: 

T v = v a v v - 5 v t 

fj, u iifa u fj,^' 

Let Y = (Y^) be a vectorfield on M and let Z = (Z a ) be a vectorfield on J\f. The Noether 
current corresponding to (Y, Z) is 

• l {Y,Z) -~ Va ^ ± v 1 ■ 
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Let J = *j be the Hodge dual of j with respect to the metric g. Noether's Theorem [37] states 
that, for any solution / of the Euler-Lagrange equations, 

d(J o Sf) = K o Sf, K:=£zL — £yL, 

where C denotes the Lie derivative operator. In particular, if L is invariant under the Lie flow 
generated by (—Y,Z) on the velocity bundle V, then Jo s/ is a closed (m — l)-form, so that its 
integral on any closed m — 1-dimensional submanifold S of M is a homological invariant i.e. depends 
only on the homology class of S. 

As particular examples, consider the case where Y is a Killing field of (A4, g), i.e. Cyg = 0. If 
the Lagrangian density I is invariant under the flow of Y it then follows that Lyh = 0, so that the 
corresponding Noether current j having Z = is divergence-free: d^j^ Y0 ^ = 0. Similarly, if Z is a 
Killing field for (Af, h) and £ is invariant under the flow of Z, then once again one gets a divergence 
free current, namely <9 M j^ z ^ = 0. 

2.2 Harmonic maps 

Definition 1. A harmonic map f is a critical point, with respect to compactly supported variations, 
of the action A, where L = I e[g] is the following Lagrangian 

:= ig^aOh^)^. 

Therefore Lo sj — |tr g /*h. 

The harmonic map action is clearly invariant under the isometrics of the domain M. and the 
target M . Thus any Killing field of either of these manifolds will yield a conservation law for the 
harmonic map. Consider in particular a Killing field Z for the target M . The corresponding Noether 
current, when evaluated on a solution section s/, is 

f = P £z a = g ^d„f b h ab z a = fsrfa = {<fly 

where <f> is the pull-back under / of the 1-form £, namely <j> = /*£, and £ = Z 9 is the dual of the 
vectorficld Z with respect to the metric h, namely — h a bZ a . 

2.3 Lie Groups 

Let G be a Lie group and U an open domain in 1", n = dimG. Given a (suitably regular) 
parametrization g : U — > G, the Lie-algebra- valued connection 1-forms w = g~ x dg and w' = —dgg~ l 
are called the Maurer-Cartan left- and right- invariant forms for G. A left-invariant form w gives 
rise to a left-invariant metric on G in the following way: Let {X a }"=i be a basis for the Lie algebra 
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g of the Lie group G. Thus w = C a X a where C a G /\ (U) arc 1-forms, with C" = Qdx' 1 , for the 
local coordinates = x G U. One computes 

itr (™ 2 ) = i^C'tr {X a X b )dx»dx v = ^Cdx^dx" =: \ v dx"dx v . (4) 

Here 77 is the Killing-Cartan quadratic form on g: n a b := ^tr (X a Xb) = hC^ d C^ c , where C^ b are the 
structure constants of the Lie algebra, defined by [X a ,X{,] = C c ah X c . Note that n is non-degencrate 
precisely when g is semisimple: in this case, one sees that the tensor h defined above provides a 
non-degenerate quadratic form on the tangent space T g G for any g G G, and thus turns G into 
a pseudo-Ricmannian manifold, of signature ((dimp)+, (dimfi)— ), where 4,p are as in the Cartan 
decomposition of g (see below) . 

2.4 Harmonic maps into Lie groups 

Consider now the case where the target manifold (TV, h) of a harmonic map / : M — > N is a Lie 
group, and h is the invariant metric defined in (4). Let W := f*w denote the pullback of the 
Maurer-Cartan form w under /. Thus W — tfXj where (p 1 — f *( ! as before, and I is a counting 
index (not a component index). Since both d and A are covariant under pullbacks, so is the equation 
dw + wAw = 0, thus dW + W A W = where d now denotes exterior differentiation on the domain 
A4. On the other hand, each dual vectorfield Zj = (C 1 )^ is easily seen to be a Killing field for 
the metric h, and by Noether's Theorem gives rise to a divergence-free current jj = (</> / ) ti . Thus 
W" = (^^Xj is also divergence free. The system of equations for a harmonic map can in this way 
be recast into the following (nonlinear) Hodge system for a Lie-algebra-valued connection 1-form 
We/\\M,Q): 

dW + W A W = 0, SW = 0, (5) 

where 5 = *d* is the divergence operator, and * denotes the Hodge dual with respect to the 
domain metric g. It is the above Hodge system that, in situations where the domain is effectively 
two-dimensional, becomes the starting point of the quest for a Lax pair, which in turn allows the 
inverse scattering method to be applied. 

2.5 Symmetric Spaces 

Let H be a complex semisimple Lie group having Lie algebra f). A real form of f) is a Lie subalgebra 
g of f) such that the complcxification of g is isomorphic to f), i.e. gc = f), so that every Xel) can 
be uniquely written as X = X\ + iX-i for some X\, X2 G g. It is always possible to realize g as the 
fixed point set of a conjugate- linear involution preserving the bracket on (): 

7-.: !>->{), tI{X) = X, T,(aX) =dT,(X), Ie|,aeC. 
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Note that for an involution r on H, the induced involution t* on f) denotes the differential at 
the identity e £ H, namely r* = dr e . In general, f) may have several non- isomorphic real forms 
arising from different choices for r. Then for G denoting the fixed point set of r in H, there is a 
corresponding Lie subalgebra of t), namely 

g = {X et)\r*X = X}. 

If g is semisimplc, then it has a maximal compact subalgebra i. This subalgebra may also be 
realized as the fixed point set of a (real linear) involutive automorphism on g. Suppose cr* is such 
an involution so that 

{ = {X e Q | a^X = X}. 

Restricting our attention to a particular real form g with maximal compact subalgebra 6, one 
may complexify t to 6c and extend <r* by complex scalars to ac, to obtain the diagram in Figure 3; 
lines between two sets indicate an involution defined on the larger set which fixes the smaller set. 



I) =flc 




Figure 3: r* and ac are in bijective correspondence 

If g is semisimple, it is known that there is a bijection between the conjugate-linear involutions r* 
and complex-linear involutions ac- Furthermore, the associated a € Aut H commutes with r (i.e. 
ar(h) = ra{h) for all h e H) [39]. For a specific example, see Figure 6, Section 4.2. 
Next, define p to be the —1 eigenspace of cr* in g, namely 

p = {X eg\a«X = -X}. 

It is known (e.g. [7]) that G/K is a symmetric space, and that g has (Cartan) decomposition 
g = fi + p, with [t, t] C 6, [t,p] C p, and [p,p] C t. Since G is semisimple, the Cartan-Killing form 

?7 (x,y) = ^c[ j x^, i,y eg 

is non-degenerate. If K is a maximal compact subgroup of G, then i] is negative definite on t and 
positive definite on p. Moreover, t and p are orthogonal subspaces with respect to 77, and { is a 
maximal compact subalgebra of g. 
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2.6 The Iwasawa Decomposition 

Using the subspaces t and p, we shall write down the Iwasawa Decomposition of g and use it to 
establish a quadratic constraint on the symmetric space G/K. 

Let a denote a maximal subspace of p that is an abelian subalgebra of g. The dimension of a is 
called the split rank of g. Let A a be the root system of the pair (g, a), i.e. 

A a = {Aea* | A ^0,0^0}, g x a ={Xeg \ [X, Y] = X(Y)X, V Y e a}, 

where a* is the dual vectorspace to a. Then A is split by the Cartan involution, i.e. A = A~ U A+ 
and the involution maps one of these sets to the other one. Introducing the nilpotent subalgebras 

n ± = Oi, 

AeA± 

the Iwasawa decomposition of the gisrj = n0a0£, for n = n~ or n + . It lifts via the exponential 
map to a decomposition for the group G — NAK, where K is the set of fixed points of a in G, and 
A, N are the subgroups obtained by exponentiating the algebras a and n - (or n + ). 

Let * denote the twisted conjugation induced by a on G, i.e. g*g' = gg'cfig) 1 and let S denote 
the orbit of the identity e under i.e. S := {ga(g)^ 1 | g G G}. Then S is a totally geodesic 
submanifold of G which is isomorphic to the symmetric space G/K under the isometric embedding 

C : G/K — > G 

gK — ► go-(g)- 1 . 

The mapping C is known as the Cartan embedding of the symmetric space into its Lie group [11]. 
Notice that if we view C as a mapping from G to G (composing g 1— > gK with C) and suppose 
q = C(g), then qa(q) — e. Thus the image of G/K in G under the embedding consists of elements 
satisfying a quadratic constraint: 

G/K = {qeG\ qa(q) = e}. (6) 

Furthermore, by the Iwasawa Decomposition, for each g e G, there exist unique k G K , a e A and 
n G N such that g = nak. Thus, C(g) can be expressed as 

C(g) — (naA;)<r((naA;)) _1 = naka(kr 1 a^ 1 n^ 1 ) = (na)ka(k)^ 1 o-(na) -1 = C(na), 

so that the mapping C is in fact well-defined and one-to-one on the solvable subgroup S of G 
consisting of elements s <G G that can be written as s = na for some a G A and n € N; that S 
is a subgroup follows from the fact that the commutator of A and N (i.e. elements of the form 
nan~ 1 a~ 1 ) lies in N. 
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Now since the two subspaces t and p of q are orthogonal with respect to the Killing-Cartan form 
r], the restriction of rj to p provides the symmetric space G/K with a natural metric. Since the 
Cartan embedding is totally geodesic, this metric agrees with the metric h induced on S = C(G) 
as a submanifold of G [28, 27, 30]. Since the Cartan embedding kills the if-factor in the Iwasawa 
decomposition, it is possible to compute this metric for S using only paramctcrizations of the 
subgroups A and N of G: Given paramctcrizations a(u) G A and n(v) € N, one computes first 
s = na and q = scr(s ), and then w = q~ 1 dq (or —dqq^ 1 as the case may be), from which h can 
be computed. 

2.7 Harmonic maps into symmetric spaces 

Because of the Cartan embedding being totally geodesic, any harmonic map into G/K is a harmonic 
map into G, and likewise any harmonic map into G whose image is contained in the submanifold 
G/K is a harmonic map into G/K. Thus the task of constructing harmonic maps into a symmetric 
space, i.e. a solution of the nonlinear sigma-model, can be simplified by reducing it to finding a 
harmonic map into the corresponding Lie group, that is to say, finding a solution of the principal 
chiral field model (see Figure 4). 




G/K 

Figure 4: Harmonic maps as chiral field and sigma models 

Furthermore, harmonic maps into symmetric spaces G/K enjoy a large group of symmetries, 
since the full group G, which could be of considerably larger dimension than G/K, acts on it 
isometrically, thereby providing a large set of conserved currents for the harmonic map. It has been 
suggested long ago [20, 51] that every field theory that can be formulated in terms of a harmonic 
mapping from an effectively two-dimensional domain manifold into a symmetric space, is completely 
integrable, and that the inverse-scattering technique can be utilized to generate new solutions from 
known ones. We now establish this conjecture for axially symmetric harmonic maps. 
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3 Integrability of axially symmetric harmonic maps 

Let G C GL(n, R) be a semisimple Lie group, and let K be a maximal compact subgroup of G. Let 
M. = R 3 have a Euclidean metric with line element ds 2 . = dp 2 + p 2 dip 2 + dz 2 given in cylindrical 
coordinates (p, z,<p) € x S 1 . Suppose / : A4 — > G/if is an axially symmetric harmonic map and 
let q be a parametrization of the symmetric space via the Cartan embedding C : G/K — > G, so that 
<j = f(p,z). By the discussion in Section 2.3, the Maurer-Cartan form w — —dqq^ 1 e /\^{G/K) 
has a corresponding pullback form W = f*w e A^-M) which satisfies the Hodge system (5). The 
divergence of this axially symmetric 1-form W in these coordinates is 

SW = - _l_ ^(^deti|g^I^) = — [d p ( P W p ) + 8 z (pW z )], 
V|detg| P 

so that the equation SW = is equivalent to d(p *2 W) = where *2 is the Hodge star on R^_, i.e. 
(*2W) P = —W z and (*2W) Z = W p . Since = and the other components of W are independent 
of <fi, we can also replace d, the exterior derivative on R 3 in these equations, by cfo the exterior 
derivative on R^_. Dropping the "2" subscripts altogether, we rewrite the Hodge system (5) as 

dW + W AW = 0, d{p*W) = 0, (7) 

where the domain is now understood to be the right-half plane in R 2 with coordinates x = (p, z). 

3.1 Lax system on a Riemann surface bundle 

The goal is to describe a Lax pair, that is to say, an over-determined linear system of equations, 
for which the integrability condition is the VK-system (7) above. We begin with the equation 
W = —dq q^ 1 and rewrite it as 

dq = -Wq, (8) 

noting that we now view q as g(x). Following [8, 5, 20], we generalize (8) to the following linear 
system 

' = -ft* 

(9) 

*Ia=o = 9 

for the unknown * : C x R^ — > C nx ™. Here, D and are generalizations of d and W, respectively: 
d 1 

D p :=d p -uj p — , uj^ := dm ^ dx d ^ fl^ := aW^ + bp^W)^, M = l,2. (10) 

The parameter A appearing in (9) is in the Riemann Sphere C, and w(X, x), a(A, x), b(X, x) are three 
C- valued functions on C x R^_. This particular form of the Lax system using three functions was 
first considered in [20]. The functions w, a and b, are assumed to be rational in A, with coefficients 
that are smooth in x, and are subject to further restrictions. Although it is possible to proceed at 
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the level of generality appearing in (10) for quite a while, in the interest of clarity we restrict our 
attention to specific choices for these functions. In particular, w, a and b shall be chosen in such a 
way that D can be viewed as a covariant derivative (on an appropriate space) and that D agrees 
with d on A = 0. 

We turn our attention first to the function w. For x = (p, z) e let 1Z X denote the Riemann 
surface which is the zero-set of the quadratic polynomial 

F X (A, w) := A 2 - 2A(z - w) - p 2 , 

i.e., 1Z X = {(A, raj) € C 2 | F x (X,w) — 0}. For p ^ 0, this is a non-singular Riemann surface, with 
branch points at w = z ± ip, away from which w is a two-to-one function of A given by 

^(A,x) = ^+z-^. (11) 

Equivalently, for a fixed zu, there are two charts on 1Z X , corresponding to the two roots X x (zu), X' x (w) 
of the quadratic in A, each defined on a slit plane C\{z + it\ — p<t<p},X x taking its values inside 
the disk |A| < p and X' x outside of it, since A X A X = —p 2 . There are two points over w — oo, and 
thus the two-point compactification 1Z X of 1Z X is the Riemann sphere C. The mapping T x : C — > C, 
T X (X) = —p 2 /X is a deck transformation on the universal cover of 7Z X and T 2 = id; since 1Z X is 
parabolic, it is known that the universal cover is the complex plane C. Note that by construction, 

w Ia=o = °- 

With this notation in place, we make some remarks concerning the various bundles which are 
involved in the study of Dty = — in (9). 

First note that the domain of the solution ^ is precisely the Riemann surface bundle B = 1Z X . 

Since the target of ^ is a Lie group H, the tangent bundle of the target has a canonical trivialization, 
which one associates with the Maurer-Cartan form taking values in rj. Using the pull-back under ^ 
of this form to B, the operator D can also be viewed as a connection on the pull-back bundle ^ _1 TiJ, 
whose fibers are isomorphic to the Lie algebra f), and the one- form O can be viewed as a section 
of this pullback bundle. The Maurer-Cartan equations in the target give rise to zero-curvature 
equations for fi in the domain. We would like to realize O as the (pulled-back) Maurer-Cartan 
form of some group element. From this perspective, we have a zero curvature result analogous to 
Theorem 2.1 of [46]: 

THEOREM 1. Let H C GL(n,C) be a complex Lie group with Lie algebra 1). LetlA be a simply 
connected domain in C x M. 2 ^, let SI = A(X,x)dp + i?(A,x)dz be a smooth t)-valued 1-form defined 
onU, and suppose D is defined as in (10). Then the equation 



13 



for ^ :U — > H has a solution iff the curvature of the connection fl vanishes, namely 

Dn + riAn = o. (12) 

Proof. Assume there exists such that Dj^> = — fij'J/, j = 1,2. Using the definitions of D and w 
in (10), it is easy to directly verify that [Di 7 D 2 ] = 0. As a result, 

= DiD 2 V - D 2 D^ = (L>2^i - £>i0 2 + 2 Oi - Oi0 2 ) 

Thus, the curvature of the connection vanishes: DO + flAfi = 0. 

Conversely, let Q(\, x) = Adp + Bdz be a connection 1-form, with A, B : U — > f), such that (12) 
holds. We first observe that ^ is to be constructed as a mapping from the Riemann surface bundle 
B into the group H. Denote by f2 the pull-back of the Maurer-Cartan form on rj to the bundle B, 
noting that it satisfies the equation df2 + f2 A il = 0. We claim that under the appropriate choice 
of coordinates on each fibre 1Z X of B, this equation is, in fact, the zero curvature equation (12). In 
particular, we see that the operator D can be expressed in terms of d, when m is chosen as the 
coordinate on 1Z X : Given \I>(A, x), define ^>{m, x) = ty(\(zu, x), x). Then 

0„*(«7, x) = ^(A(tu, x), x)<9 M A + a„* - *( A ( ro , x), x) + 0„* = D„*(A, x). 

OA OaW CM 

Consequently, the zero curvature condition (12) can be rewritten for :— Cl(X, x) as 

g?0 + flAO = 0. 

At this point, we may appeal to the standard zero-curvature theorem for f] (see, for instance, Cor. 
4.24, p. 81 in [25]), to conclude there exists a mapping F :U — > H such that f2 = —dFF^ 1 (and 
consequently = —DFF^ 1 ) locally. Since U is simply connected, this statement holds globally. 
Calling the global mapping ^, the converse direction is proved. □ 

We note that the mapping F above will not be unique, since one may equally well consider 
its right-translation by a constant group element h e H: (D(Fh))(Fh)^ 1 = (DF)hh~ 1 F~ 1 = 
(DF)F' 1 . 

Finally, we would like to recover dq = —Wq in (8) from (9) by further restricting the functions 
a and b such that O| A=0 = W. In conjunction with the zero curvature condition, the system of 
equations for a, b in (10) are thus found to be [20]: 

«Ia=o = !> & Ia=o=°> a 2 +p 2 b 2 -a = 0, Da = p*Db. (13) 
For subsequent sections, it will be useful to fix a and b. We do so now, in the form of a Lemma. 
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Lemma 1. Let fi p := aW^ + bpfoW)^, fi = 1, 2 and define a, 6 to &e 

p 2 A 
a ( A ' x ) = ^2^2> 6 ( A ' X ) = ]^^2- 

Suppose W a 1-form on M.\, then the zero- curvature condition (12) is satisfied if and only if the 
Hodge system (7) is satisfied. 

Proof. One may check directly that a, b defined in the Lemma satisfy the necessary requirements 
stated in (13). Expanding (12) using the definition of fi, we obtain 

a(dW + W A W) + b d(p * W) = 0, 

since W is independent of A. For fixed x, this equation must hold for all A. Since a and b are 
linearly independent functions of A, both equations in (7) must hold simultaneously. □ 

3.2 Gauge freedom 

By the general definitions of lo,zu in (10), it is easy to see that the covariant derivative D has the 
property that it vanishes on zu, and consequently on any sufficiently smooth function of zu: 

Proposition 1. The kernel of the operator D consists of arbitrary matrix-valued C 1 functions of 
U7(A,X) = £+z-±. 

Proof. This follows from Implicit Function Theorem. Observe that F(p, z, A, w) := p 2 + 2X(z — 
w) — A 2 = and thus away from the branch points w = z ± ip, where d\F = 0, we have that 
A = \(p,z,m) is C 1 and d p X = -g^ = -£g, and similarly d z X = -^f . Let / : C x -> C 
be a C 1 function with Df = 0. Then g{p,z,w) := f(p,z,X(p,z,w)) is also C 1 , by chain rule 
dp9 = 9 p f + d\fdpX = and similarly d z g = 0. Thus g is only a function of zu. Applying this to 
each entry of a matrix C € ker D establishes the result. □ 

Thus, if <]/ solves (9), then so does 

*'(A,x)=$(A,x)CMA,x)), (14) 
where C : C — > G is any matrix- valued complex curve such that lim C(w{X)) = C(oo) = /. Visibly, 

A->0 

the converse is also true. 

Corollary 1. Suppose : C x M 2 ^ — > C" xrl satisfy the same linear equation 

L>* = -fi*, L>*' = -fi*'. 
Then there exists a matrix-valued complex curve C : C — >• <C nxn such that (14) holds. 
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Proof. Let V = f^ 1 *'. Then DV = -^D^' 1 ^' + = - *" 1 0'J' / = 0. By 

Proposition 1, there exists a C 1 matrix function C{w) such that \&' = *C. □ 

This is the so-called gauge freedom in the initial value problem (9). We denote by [*] the 
equivalence class of * under gauge transformations, so that e [*] iff there exists a map G : C — > G 
with C(oo) = / such that (14) holds. Note that the results of Proposition 1 and Corollary 1 arc 
not exclusive to our (fixed) choice of w. 

3.3 Reality conditions 

As seen in the above, the Lax system (9) has a plenitude of solutions, not all of which may be 
of interest to us or indeed useful for the purpose of ISM. One may ask for example if there are 
solutions ^(A) that remain in the real group G for A ^ 0. To address this question, it is necessary 
to define (following [45, 44]) the concept of G-reality. 

Definition 2. Let T> be a domain in C containing the origin that is invariant under complex 
conjugation, i.e. T> — T>. Let G be a real form of the complex Lie group H, consisting of elements 
in H that are fixed by the involutive automorphism t : H — > H . A mapping g : T> — > H is said to 
satisfy the G-reality condition if 

r{g(X)) = g(X) for all XeV. 

Similarly, a mapping £ : V — >• f) into the Lie algebra of H is said to satisfy the G-reality condition 
ifr*(£(A)) = £(A). 

With this definition we now have: 

Proposition 2. There exists a domain V as in Definition 2 such that for every solution of (9) 
there is a complex curve C : C — > C" x " for which 

r(tf (A)) = *(A)G(ro), /or a// XeV. (15) 

Furthermore, C satisfies the quadratic reality constraint 

C{w)t{C{tu)) = I. (16) 

Proof. One notes that for any fixed x, fl = aW + pb * W is a Lie-algebra-valued mapping A 
fi(A) € rj (recall e q). The functions 07(A), a(A) and 6(A) are equivariant under conjugation, 
namely w{\) = w(X) and similarly for a, b. Since r* is conjugate-linear, this implies that f2 satisfies 
the G-reality condition 

r*(Q(A)) = 0(A). 
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From here it is easy to see that D(r(\J'(A))^' _1 (A)) = 0, noting that the conjugate inside the 
expression is evaluated after differentiation. By Proposition 1, we conclude that there exists C(w) 
such that (15) holds. Taking limits as A — > 0, we observe that 

lim r(*(A)) = r(q) = q = lim *(A)C(o7) = gC(oo), 

A— >0 A— >0 

and therefore C(oo) = I. This means that there exists a neighborhood of A = for which C is 
invertiblc. On this neighborhood, apply the involution r to the equation and substitute the original 
expression for r(*) to obtain 

tf (A) = t(*(A))t(C(h7)) = *(A)C(w)r(C(ro)), 

from which we can immediately deduce the quadratic reality constraints C(w)t(C(o7)) =7. □ 

3.4 Involutive symmetry 

We now move on to the consequences of the initial data q of (9) belonging to G/K G. Since 
elements of the symmetric space satisfy the quadratic constraint (6), we have a corresponding 
symmetry for W = —dqq^ 1 and hence also for Q: 

Wq + qa*(W) = 0, fl(X)q + qa, (0(A)) = 0. 

The above now imply a certain involutive symmetry for * under the inversion T x . 

Proposition 3. There exists a domain T> as in Definition 2 such that for every solution ^ of (9) 
there is a complex curve J : C — > C" x ™ for which 

*(T X (A), x) = g(x)tr(tf (A, x)) J(w(X, x)), /or X eT>. (17) 

Furthermore, J satisfies the quadratic symmetry constraint 

a{J{m))J{w)=I. (18) 

Proof. First, observe that the functions a and 6 chosen in Lemma 1 are equivariant under T x , i.e. 

a(T x (A),x) = l-a(A,x), &(T x (A),x) = -6(A,x), 

as a result of which, O is equivariant under T x as well 

0(T x (A,x),x) = W-fi(A,x). 

Let us then define >3/(A,x) := *(r x (A),x). It is easy to check that D commutes with T x , i.e. 
£>^(A,x) = (D4 , )(T X (A), x), since vj is invariant under T x . It then follows that ^ satisfies the 
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equation = {—W + £l)*5f, and using the above mentioned symmetries of fi and W , we observe 
that the function \&' := qcr^) satisfies the same equation as ^, namely D^>' = — fl^f' (although not 
necessarily with the same initial value as Thus by Corollary 1 we must have a complex curve 
J for which = ^>'J(w). In other words there exists a mapping J = J(zu) such that (suppressing 
dependence on x) 

*(A) = qa(^(T x (\)))J(m(Xj), (19) 

for all A e C. Replacing A with T X (A) in the above yields (recall that w is invariant under the 
inversion T x ) 

*(T X (A)) = g*(*(A)) J(w{\)). (20) 

Substituting (20) into (19), we conclude that J(zu) must satisfy the quadratic constraint, and have 
proved the proposition. □ 

Note that J(vj) encodes the asymptotic behavior of ^ as A — > oo: Since T" X (A) = — p 2 /A, we 
may take the limit A — > in (20) to conclude 

lim *(A,x) = qa(q)J(oo) = J(oo), (21) 

A— »oo 

which also shows that <J/ at A = oo is a constant map. 

3.5 Inverse scattering and vesture for harmonic maps 

Suppose that go : M G/K is a given axially symmetric harmonic map into the symmetric space 
G/K, q = qo(p, z). Assume be any solution to (9) with initial data qo, let Wq = —dq a q^ 1 and 
£! = aWo + bp * W - Vesture (dressing) refers to the possibility that other solutions q — f(p, z) 
to the harmonic map may be generated by finding a matrix \ = x{\ x ) with the property that \I>, 
defined by 

*-X*o, (22) 

solves the PDE in (9). If in addition we require ^ and >1> to have the same asymptotic behavior 
as A — > oo, it would follow that 

X(oo,x)=7, Vxel^. (23) 
If such a matrix \ can be found, then setting A = 0, one obtains 

«(x)= X (0,x)gb(x), (24) 

which gives us a new solution q(x) to the harmonic map problem, obtained by "dressing" the seed 
solution <7o(x). The process of arriving at q from q via \ is thus formally similar to the inverse 
scattering discussed in the Introduction. Compare Figures 1,2, and 5. 
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Figure 5: Vesture for Harmonic Maps 

We shall now outline the procedure for finding the dressing matrix X an d generating new har- 
monic maps q via ^F. Choose an "initial seed" map go, an axially symmetric harmonic map of K 3 
into G/K, and solve the linear system (9) to obtain an initial matrix solution of the Lax system, 
*o- Let 

O := aWo + bp * W 

where Wo := —dq^ 1 . Now given any invertible matrix- valued map X = X(A, x) which is rational 
in A and smooth in x, we may define the one-form by 

O:=(D X - X O ) X - 1 . (25) 

By construction, then, 

D X = X O - X . (26) 

We are of course only interested in those X for which g(x) = \l/(0,x), the result of the dressing 
(24), still belongs to the symmetric space G/K. In particular, X needs to satisfy symmetries 
corresponding to those of <J/, e.g. det'I'(A) = 1 for all A e T>. On the other hand, using the 
definition of O, 

D{det X ) = det X tr(x _1 Dx) = -det X tr( X _1 [xO - O x }) = -det X trfi. 

A priori, tr O does not vanish, so we cannot conclude that det X = 1. To address this, we make the 
following modifications (this is analogous to the modifications in [8, 9].) Set 

n' = fi - itr fi I nxn , X ' = (det X )^ 1/n X, (27) 
n 

so that trO' = and det X ' = 1- Observe that the following still holds: 

D X ' = x 'O - O'x'. (28) 

This means that VP' = x'^o solves (9) with 0' in place of O. Since tr O' = implies that O' e f), 
we now have that f'eff and thus g(x) = *'(0, x) e i?. 

In order to further ensure that q & G/K ^ G and therefore it is indeed a new map into the 
original symmetric space, further restrictions need to be imposed on X , namely the G-reality and 



*U=o =<?o 



A=0 



^*o(9o,A) 

*=X*o 

-*(«,A) 
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involutive symmetry conditions corresponding to (15) and (17). In the first case, r(x(A)) = x(A). 
For the second condition, we simply observe that if x(A) is a solution to (23)-(26), so will be 
X'(A) = qa(x(T x (X)))a(qo) , where q = x(0)go- We are thus going to require that x' = X t° ensure 
the resulting ^ possesses involutive symmetry as well. 

The above restrictions will guarantee the map q is in the symmetric space G/K. In order for q 
to be a harmonic map, however, still further restrictions arc needed, this time on the pole-structure 
of x as a rational function of A. To summarize, 

Definition 3. A mapping x : C x — » GL(n, C) is a dressing matrix for q , an axially symmetric 
harmonic map of R 3 into G/K, if it satisfies all of the following: 

X(oo) = I (29) 
r(x(A)) = X(A) (30) 
X(A) = qa( X (T x (\)))a(q ), for q := X (0)q . (31) 

In addition, the poles of \ are to be restricted in such a way that the poles of £1 defined by (25) 
are precisely those of £1 , and with the same residue at each pole. 

We now have the following theorem. 

THEOREM 2. If qo is an axially symmetric harmonic map o/R 3 into G/K and \ is a dressing 
matrix for qo, then q = x(0)<7o is also an axially symmetric harmonic map o/R 3 into G/K. 

Proof. Given qo,x, construct O as in (25). Modify O to 0' as above, noting that this modification 
cannot introduce any new poles. Given a, b defined in Lemma 1, set 

W := ft' - lp(*tf). 

With this definition of W, one has 

n' = aW + bp * W. 

We now claim that, for fixed x <E R^, W(-,x) is holomorphic on the Riemann sphere, and is 
therefore constant in A, by Liouville's Theorem. Recall that Slo is holomorphic on the Riemann 
sphere except for simple poles at A = ±ip. From the definition of W it is clear that the only possible 
poles for W can be at A = ±ip. Recall also that the residue of the meromorphic one-form a on the 
Riemann surface R at a point p is defined as 

ReSpa = Res z=p / = J a, 

where a = fdz in local coordinates about p, and 7 is a small circle around p. Calculating the 
residue of W at its only two poles, A = ±ip, we obtain 

Res x = ip W = Res x =iptt' - i(Res x=ip * O') = Res x=ip fl - i(Rcs x=ip * ft), 
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since \p is holomorphic at A = ±ip. Next, from the definition of \ we know that the poles of Slo 
and agree (with the same residues), thus 

ResA=i P fi = Res\ =i ptt = W Res\= ip a + p * W Res x =i P b = -^W + ^ * W . 



The final two equalities follow from the fact that Wo is independent of A as well as using the 
definitions of a and b, respectively. Combining this result with the definition of SI', we have 



i.e. the residue of W at A = ip vanishes. Repeating the calculation now for A = — ip we conclude 
that W has no poles anywhere on the Riemann sphere, which establishes the claim. 

We thus have O' = aW + bp(*W), where W is now independent of A. Since ^ = x'^o is a 
solution to the Lax system (9) by construction with fl' in place of SI, we may apply Theorem 1 to 
conclude that SI' has zero curvature. By Lemma 1, this implies W satisfies the Hodge system (7). 
Thus q is indeed a harmonic map. The symmetries of \ 011 the other hand ensure that q € G/K, 
and we are done. □ 

3.6 Algebraic solution of the vesture problem 

In general, the task of finding a dressing matrix \ for a given q is seen to be equivalent to a 
Riemann-Hilbert problem in the complex plane [8] , the complete solution of which requires solving 
a certain system of integral equations coupled to algebraic equations. In some cases, however, 
referred to in the literature as "the solitonic case" , a special Ansatz for \ can be employed which 
guarantees that the Riemann-Hilbert problem has a trivial solution, and the aforementioned system 
then reduces to a purely algebraic system of equations. Such an Ansatz was considered in [8, 9, 20], 
etc., postulating x to be a rational function of A with a number of prescribed simple poles: 



We say an n x n matrix A has a simple pole at A = Xq if (A — Xq)A is non-zero and analytic 
at A = Aq. The poles Afe = Afc(x) are prescribed according to a certain recipe given below, and 
Rk = -Rfc(x) are matrices to be found subsequently. This Ansatz reduces the linear system (26) to 
a set of algebraic equations for a collection of singular matrices Mk, to be defined shortly. We are 
assuming that a solution ^o(A) of the Lax system (9) corresponding to a harmonic map qo is given, 
and that the poles Afe being prescribed are not already poles of ^o- Therefore the new solution ^ 
constructed by way of the dressing matrix \ will have poles at Afe. The example we will consider 
later shows that this procedure will lead to an axisymmetric harmonic map with ring singularities. 



ResA=i P *Sl = — * Wo — —Wq = —iRcs\ = i P Sl, 




(32) 
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We also observe that if x has a pole at A = Afe, then by (31), x must also have a pole at 
A = Tx(Afc), while by (30), must have a pole at A = Afe. The above mentioned symmetries also 
suggest that poles Afe(x) are in fact being prescribed on the Riemann surface 7?. x , in the following 
sense: Let {^k)k=i ^ e ^ distinct non-real complex numbers 2 ; without loss of generality, we may 
choose these values to be in the upper half plane. Define the poles 

Afe = A x (tn fe ), Ajv + fc = X' x (wk), k = 1,...,N, (33) 

where A x , A x are the two charts for 1Z X . In particular, the poles come in pairs related by the deck 
transformation on TZ X : 

Aw+fe = T x (A fe ). 

Moreover, the definition of A fe in terms of the two charts A x and A x implies that raj(Afc(x), x) = ju k , 
which upon differentiation with respect to x implies ^j-d x Afe + rf x cc = 0, in other words, 

rfAfe +u;(A fe (x),x) =0. (34) 

Now, since XX 1 = I holds for all A, it also holds at A = Afe, which is a pole of X- It follows 
that X~H^) must be holomorphic at A = Afe and non-zero, and further, that it is singular since its 
determinant has to vanish as A — > Afe. Thus, 

RkX~ 1 (*k)=0, (35) 

and in particular, the matrices Rk must also be singular (deti?fe = 0). We may use the symmetry 
properties of \ to find what they imply for the Rk . For the sake of brevity, let us denote 

tf^(x) := lim ^(A.x), 

A^A fc (x) 

and similarly for V^o/c and Xfc^ 1 - 

First we observe that the following symmetry of x is implied by (30) and (31) (again, suppressing 
the x dependence): 

X- 1 (A)g = g <7(r[ X - 1 (T x (A))]). (36) 

Multiplying on the left by Rk and taking the limit as A — > Afe , we use the numbering convention in 
(33) to obtain 

= i?fe 9 oTr[x- 1 (r x (A))] A=Afc = J? fe go<7T[x- 1 (A)] A=Ajv+fc , (37) 

2 Our setup does not allow for these poles to be real. We should mention however, that there is a closely related 
setup which does allow the poles to be real, generalizing the approach in [8] . This is in particular important in the 
context of the Einstein Equations, if one insists on obtaining solutions which have black holes, as opposed to naked 
singularities. We will pursue this line of inquiry in a future paper. 
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where the addition in N + k is mod 27V, i.e. for k > TV ', N + k = k — TV . Another symmetry of \i 
from (30), is 

X - 1 (A)=r( X - 1 (A)), (38) 
which upon left multiplication by R k and taking the limit yields 

= R k r[ X - 1 (\)] x=Xk . (39) 

Equations (37)-(39) appear to be an over-determined system, since there are two matrix equations 
for each Rk- In order to be able to reduce the number of equations, we now make an assumption 
about the involutions r, a, namely that they should be given by conjugation with respect to the 
same element. 

Assumption 1. Assume that there exists an element r e H with T 2 = /, such that 

r{g)=Y{g*)- 1 T, v{g) = YgY. (40) 

We remark that although there are semisimple Lie groups possessing involutions which cannot 
be realized in this way, the assumption can also be rephrased in terms of inner equivalence classes in 
order to identify those Lie groups which are allowable [4, 39]. In particular, any Lie algebra having 
a Dynkin diagram with no symmetries satisfies this restriction. If on the other hand, the Dynkin 
diagram possesses nontrivial symmetries, such is the case for example if H = SL(n, C), then there 
will be real forms of H where the corresponding involution is not realizable as conjugation with 
respect to an element, e.g. G — SL(n,M.). However, even in those cases, other real forms of H 
may still satisfy Assumption 1. As an example, we note that such a matrix Y exists for the pseudo- 
unitary groups SU(p, q), to be discussed in Section 4.2. Under this assumption, the overdetermined 
system (37)- (39) can be re-expressed as 

Rkqo [x&N+k)]* = (41) 
R k Y[x(X k )]*Y = 0. (42) 

Or, using the Ansatz in (32), 

2N 

R k q + J2\ ^RkqoR* - (43) 

j=1 A^v+fe - Aj 

2N 

Rk+J2\ T RkYR*Y = 0, (44) 

j=i Afe - Aj 

where the addition in N + k is mod 27V, i.e. for k > N, N + k = k — N. 

Now clearly if RN+k = HkRkQoY for some invertible matrices H k , then the two equations above 
would become equivalent. We will use Proposition 3 to show that it is always consistent to make 
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such a choice. For k = 1, . . . , 27V, define the matrix 

M k := ifc(x)* fe(x), (45) 

noting that by assumption, is not already a pole of ^q. Recall that the matrices R k s are known 
to be singular, from which it follows that the M k are also singular matrices in C" x ™ for each k. 
Indeed, by (35), 

M^ 1 - i? fe *ofe*fe 1 = RkXu 1 = 0. (46) 
The symmetries of VP, (17) in particular, together with (6) and (18) imply 

tf-^x) = ( r[J(tn(A))] ( r[vI/- 1 (r x (A),x)] ( 7( g ). 

Multiplying by M k and taking the limit A — > A^ we obtain 

M k a{J{w k ))a{^-\ k )c7{q) =0 => M k T J{w k )^ N \ k T = 0. 

On the other hand, by (46) we also know that 

M N+ k^- N \ k = 0. 

Comparing the two equations above, matrices M N+k and M k T J(w k ) have the same null space, and 
thus it is consistent to assume that 

Mjv+fe = H k M k TJ{vo k ) = H k M k a(J(m k ))T (47) 

for some invertible matrices H k . Note that as A — > X k , the complex curve J(w) appearing in 
Proposition 3 also satisfies J '(w n +k)J '{^k) = Inxn- Consequently, this yields the relation between 
i?jv+fe and R k , making the two equations in (41) coincide, since the symmetry of in (17) implies 

R N+k = M N+k ^ N+k = H k M k YJ{w k )%^ N+k = H k M k %^q T = H k R k q T, 

as claimed. It is thus enough to keep only the equation (44), and rewriting it in terms of M k we 
arrive at the following (nonlinear) system 

27V 

M k% k F + J2\ T M k S kj M*=Q, k = l,...,2N, (48) 

j=i A fe - *i 

where we have set 

s kj := * - fe 1 r(^)- 1 = r*5 fc r* 0j r. (49) 

It remains to ensure that the matrix \ thus defined has the correct pole structure as in Definition 3. 
In order to do so, we simplify (48) further using an additional reduction: We ssume that the R k 
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(and hence the M k ) are rank-one matrices, i.e. there exists non-zero vector functions Ufe,Vfe for 
fc = l,...,2N such that 

M k = u k v* k . 

This rank-one assumption is consistent with detx having only simple poles at the X k , and is 
equivalent to assuming that the matrices Xfe 1 have rank n — 1. 
From (47) we observe that it would then be consistent to take 

Vfe = TvAr+fc, k=l,...,N. 

Thus a complete set of unknowns for the problem are the vectors {ufc} 2 .^ together with {v k } k=1 , 
i.e. 3iV vectors in C", instead of the 2N matrices M k in C" x ". Together they satisfy 

2N 

E 1 T <S kj v jU * = -v^o^r (50) 

Moreover, as a consequence of Theorem 2, we have the following 

THEOREM 3. The dressing matrix \ can be constructed from (50) in such a way that 

has the exact same pole structure as S7o • Furthermore, without loss of generality, vectors {vk}^! =1 
can be taken to be arbitrary constants, i.e. independent o/x. 

Proof. Clearly, it is necessary for f2 to be holomorphic at A = Afe for any k, and any x^xt, where 
the points Xfe are where Afe = ±ip (the poles of a and 6, and thus of Clo). It is easy to see that for 
any choice of zu k = a + i(3, setting z — a, p — will yield such a point Xfe. Thus the analysis 
below can be carried out at all but finitely many points in R+. 

We use the Ansatz (32) to rewrite the right-hand side of (25) in terms of inverse powers of A — Afe 
and perform a residue analysis. It turns out, thanks to (34), that the coefficient of (A — Afe) -2 in 
fact vanishes at A = Afe, and we end up with 

2N 

-Cto + E \ _ \ (dRk — Rk^o) 



{D X -X^)X^\ X= 



A=A fc 



fc=l A " Afe 



X- X (Afe). 

A— X k 



Since fio is holomorphic at A = Afe (and away from x = Xfe), we obtain that the coefficient of 
(A — Afe) -1 also has to vanish at A = Afe. Thus 

[(iR fc (x) - i?fe(x)O (Afe(x),x)] X " 1 (Afe(x),x) = 0. (51) 

Suppressing the x-dependence for brevity, we have 

dM k = dR k ^ ak + i?fc(L»* )(Afe) = [dR k - R k n (X k )}^ ok 



25 



and using the rank one assumption, (51) becomes 

= dM^Xk 1 = dMk*- 1 = du^- 1 + Ufcdv£ttfc K 

Now by (46), M^^ 1 = implies that v^^ 1 = 0, which in turn gives dv^^ 1 = from the above 
equation. Since VE^ 1 = 'J'ofc X& 1 is rank n — 1, dvfe has to be a multiple of Vfe, i.e. one must have 

dv fe = ftvfe, (52) 

for some scaZor-valued 1-form ft. Applying d now, we obtain 

= dhw* k + ft A ftv£ = dftv^, 

which, since are nonzero, implies dh — 0, and hence by Poincare's lemma, ft = dj for some 
function 7. Thus the differential equation (52) can be solved to get v£ = e 7 v£(x ), i.e., the vector 
functions are scalar function multiples of a fixed vector Vfco- From (50) it follows that if 
are scaled by a factor of e 1 , then Mfc and the new solution q constructed by this method will be 
independent of the function 7, and therefore without loss of generality one can assume that v^'s 
are constant to begin with. □ 

The equations for Mk (48) now become a linear system for the unknown vector functions U& in 
terms of the constant vectors v&, using (49): 

°*J' U 3 = b fc' ak J := . 1 = v^5fcjVj, b* k := -vjtfj^r. 

These equations can be written in matrix form as 

AU*=B*, (53) 

where U = U(x) is the (n x 2N) matrix whose columns are the vector functions , A = A (x) = (a^ ) 
is a 2N x 2N matrix function, and B = B(x) is the n x 2N matrix whose columns are the vector 
functions defined above. 

If the matrix A can be shown to be invertible 3 , at least in a neighborhood in then the 
above system has a unique solution U* = A~ 1 B* in that neighborhood. From there, one can then 
calculate the matrices Rk and the dressing matrix x, and setting A = 0, the new solution q is found 
to be 

2N 

<z(x) - <zo(x) - J2 T-TT Ufe W v fc*ofc( x )9o(x). (54) 
fc=i AfeW 

We will see in the examples of the next section, that the matrix A in general is not invertible 
everywhere in the domain K3_, and indeed the zero set of detA has a geometric significance for 



3 As far as we know, this point has not been addressed in previous studies of the dressing technique. 



26 



the dressed harmonic map o/(x). However, we can prove that under appropriate conditions on the 
arbitrary vectors Vfc, the matrix A(pc) is in general invertible for large |x|, i.e. in a neighborhood 
of infinity. First we recall the following: 

Definition 4. A matrix A £ C nx ™ is called strictly diagonally dominant if 

J~] \a,ij\ < \au\, for alii = 1, ...,ra. 

By the Levy-Dcsplanqucs Theorem (see e.g. [43]), the determinant of a strictly diagonally 
dominant matrix is non-zero. We prove 

THEOREM 4. Let vi, . . . vat eC" be N arbitrary complex vectors satisfying the condition that, 
for all k = 1,...,N, 

EKM<^K r ^l- ( 55 ) 

Set VN+j = Tvj /or j = 1, . . . , N. Then, there exists R > smc/i that for all x £ wii/i |x| > i?, 
i/ie matrix A(x) £ C 2Afx2Ar wii/i elements 

w/iere i/ie Aj 's are as m (33) and Sij as in (49), is strictly diagonally dominant, and hence invertible. 
Proof. Let zu\, . . . , m n £ C be TV distinct points in the upper- half complex plane, lm(wj) > 0. Set 

Wj = zj + isj. 

For j = 1, . . . , N let (rj,6j) denote TV systems of elliptical coordinates, with real parameters (zj,Sj), 
on K 2 ^. In terms of the Cartesian coordinates (p, z) we have 



p = \J r j + s] sin @j i z = z j + r j cos #j ■ 

Recall that Aj and Ajv+j are the two roots of the quadratic polynomial 

A 2 - 2(z - Wj)\ - p 2 = 0. 

In terms of the elliptical coordinates, we have 

Xj = ( rj - is^icosBj + 1), X N+j = (r 3 + ^(cosf^ - 1). 

Note that as |x| — »■ oo, we have Tj —¥ oo for all j, and thus \Xj\ — > oo as well. 
Fix k £ {1,...,AT}. Then 

|Afe - Afe| = 2s fe (l + cos6> fc ), (56) 
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while, for 1 < j < N and j ^ k, 

\X k - Xj\ = \r k (l + costffe) - r,(l + cosBj) - i (s k (l + cos9 k ) + Sj(l + cos0j))| > s k (l + cos9 k ). 

(57) 

For N + 1 < j < 2N on the other hand, 

| A* - Xj\ = \r k (l + cos9 k ) + rj{l - cosOj) - i (s k (l + cos9 k ) + Sj(l - cos#j))| > r k (l + cos9 k ). 

(58) 

It is clear that the same relations also hold for k € N + 1, . . . , 2N, except that (57) now holds for 

N + 1 < j < 2N and (58) holds for 1 < j < N. 

Now recall that 5'ofe( x ) = ^o(Afc(x), x). We have shown (Prop. 3) that, as A — > oo, the matrix 
approaches a constant matrix J that satisfies the quadratic constraint, see (21). It thus follows 

from the definition of matrices SV, that Sij — > T as |x| — > oo for all i,j = l,..., 2N, and hence 

\r k S k j\rj — > v k Fvj, as |x| — > oo. (59) 
For vectors Vj satisfying the condition (55), and for r k large enough, we have 

N , ^ , N 



Krv,l , x-vjv,- ^ |v*r Vfc | 

Therefore, as |x| — > oo, using (56)-(58), 



l v fc r ^jl < l v fe rv fcl 



Combining this with (59), we arrive at the desired result, namely the diagonal dominance of A. □ 
We have thus shown 

THEOREM 5. Let G be a real semisimple Lie group admitting a pair of commuting involutions 
which satisfy (40) and suppose K is a maximal compact subgroup of G. Then the Hodge system 
(7) for an axially symmetric harmonic map from M 3 into the Riemannian symmetric space G / K is 
integrable by way of ISM. In particular, it is always possible to generate new harmonic maps from 
any given one using the dressing technique. 

We note in passing that even though it appears that one needs to know the whole solution 
*o(A, x) to the Lax system (9) with initial data qo in order to find the matrices M k , the only 
information one needs about V&o is its value at the new poles X = X k (x), i.e. the matrices \&ofc( x )- 
These can be found by "integration" (more precisely, exponentiation) from the seed solution qo(x), 
in the following manner: 
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It is straightforward to check that for any given map 0(A,x), the map VK X ) — <K^fc( x )i x ) 
satishes (D</>) (A& (x) , x) = di/j(x). In particular, the matrices "J/o/c satisfy the following equation 

d^ok = -^ fc*ofc = -(a k W a + pb k * Wo)* fe, 

where the k index denotes evaluation at A = Afe(x) as before, and Wo we recall is the Maurer- 
Cartan form corresponding to the seed solution go- Since is an element of the group G, the 
Maurer-Cartan form corresponding to ^ok, Wk := —d^ok^ok = a kWo + pbk * Wo, is known once a 
seed solution and a set of poles Wk are specified. Using the parametrization we have for the group, 
it is thus possible to recover \J/ fc from the Lie algebra element Wk by exponentiation. 

4 Applications: Generating Solutions 

4.1 Agreement with previous results 

The techniques of this paper may be used to confirm calculations for stationary axisymmetric solu- 
tions to the Einstein Vacuum and Einstein-Maxwell equations. We briefly indicate the approaches 
taken in these two particular cases of interest in gravitation before exhibiting a general approach 
which subsumes both. The simple pole ansatz (32) for \ appearing in this paper was introduced by 
[8, 9], where stationary axisymmetry solutions to the Einstein vacuum equations were first studied. 
The authors consider a two-dimensional representation of SL(2, M), calculating AT-soliton solutions 
explicitly. Using the Minkowski metric as an initial seed, the authors recover the Kerr-NUT metric 
in the 1-soliton case. In contrast with the approach of this paper, the poles A& appearing in \ are 
chosen in conjugate pairs and thus a different constant vjk is needed for each pole A^ . 

Generalizing the approaches used in [32, 34] for the vacuum equations, the authors in [20] for- 
mulate the Einstein-Maxwell equations as a harmonic map into a symmetric space having isometry 
group SU(2, 1). The choice of a 3-dimensional representation of this group realises the symmetric 
space as the complex hyperbolic plane He = SU(2, 1)/S(U{2) x U(l)). Although explicit exam- 
ples of dressing do not appear in [20] , their discussion motivates our generalization to the complex 
Grassmann manifolds SU(p,q)/S(U(p) x U(q)). 

4.2 New results: Noncompact Grassmann Manifolds 

We shall describe in detail a solution-generation method for the harmonic map equation, in the 
particular case where the domain is R 3 and the target is a non-compact Grassmann manifold G P , q , 
a 2pg-dimensional Riemannian symmetric space realized as a quotient of the group G — SU(p, q) 
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by its maximal compact subgroup K — S(U(p) x U(q)): 

G p , q :=SU(p,q)/S(U(p) x U(q)). 

We restrict our attention further to axially symmetric solutions, the domain thus becoming effec- 
tively two-dimensional. The results mentioned above on the Einstein vacuum and Einstein-Maxwell 
equations are thus about axisymmetric harmonic maps into and £2,1, respectively. As a further 
example, the study of the particular group SU(2,2), and the associated Grassmann manifold ^2,2 
is of physical interest, as this group is the universal (double) cover of SO(A, 2), the conformal group 
of the Minkowski spacetime M 3 ' 1 [17, 47]. 

We begin by recalling (e.g. [30]) that the real semi-simple Lie group SU(p,q) can be identi- 
fied with the subgroup of SL(n, C) consisting of matrices in C nx ™ preserving a pseudo-Hermitian 
quadratic form having p plus signs and q minus signs, where p + q = n: 

G = SU(p, q) = {9 € SL(p + q, CJIIV)-^ - g}. 

Here, T — T pq is the n x n block-diagonal matrix diag(/ pxp , — I q xq)- Note that one may equally 
well choose any Hcrmitian matrix T conjugate to T p q to define the unitarily equivalent presentation 
of the group SU(p, q), so for ease of notation, we suppress the subscripts p, q on the matrix T. The 
corresponding Lie algebra to the group is 

su(p, q ) = {x e gl{p + q,C)\ + tx = o, tr x = o}. 

Define the following two involutions on H = SL(p + q, C) 

rfo) = r( 5 *)- 1 I\ a(g) = TgT. (60) 

It is easy to see that r, a are involutive automorphisms, SU(p, q) = {g € SL(n, C) | r(g) = g}, and 
K = {g G G I cr{g) = g} is isomorphic to S(U(p) x U(q)), which is a maximal compact subgroup of 
G, making G/K a symmetric space. The induced Lie algebra involutions are given by 

T,X = -TX*T, cj*X = TXT, Ie 8 . (61) 

The involution 0% complexifies to the complex-linear map ac(X) = TXT so that one has a diagram 
in Figure 6 analogous to that of Figure 3. 

The induced involutions give rise to the Cartan Decomposition g — t © p. Here, p = {X e 
g I TXT = —X}, and we observe that 6, the Lie algebra of K, is a maximal compact subalgebra 
of 0, which is indeed isomorphic to s(u(p) x u(q)). The quadratic constraint on the image of G/K 
under the Cartan embedding is 

G/K = {q G G I qTqT = I}. (62) 
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g c =sl(n,C) 




t = s(u(p) x u{q)) 



Figure 6: t*(X) = — TX*T (C-linear ) and <Jc{X) — TXT (C-linear) are in bijective correspondence 

It is straightforward to confirm that the corresponding Lie algebra, su(p, q) = g has dimension 
(p + q) 2 — 1, while dimfi = p 2 + q 2 — 1, and thus dimQ p ^ q = dimp = 2pq; further details in this case 
are carried out in [7, p. 39]. 

One also notes that since G = NAK by the Iwasawa Decomposition, each g £ G can be 
expressed in the form g = nak for unique k e K , a e A, n e N; consequently, for G = SU(p, q) 
with the above defined involutions 

gcr(g)~ = naa(a~ 1 )(j(n~ 1 ) = naTa^ 1 ^^. 

On the other hand, by the definition of G and the fact that r(g) = r(.g*) _1 r, 

gcr(g)^ 1 = gTg~^^T = gg* = nakk*a*n* = naa*n* , 

since k e K is a unitary matrix. Thus in particular the image of the symmetric space under 
the Cartan embedding consists of Hermitian matrices. Moreover, having parametrizations for the 
subgroups A and N suffices for obtaining a parametrization of the symmetric space. This will be 
carried out explicitly for SU{2, 1) in the next section. 

A crucial fact about Grassmann manifolds is that the natural embedding G p > yq > =-> G P ,q, for 
p' < Pi q' < is totally geodesic. This is due to the following more general result (see [27, Thm. 
IV.7.2]): 

THEOREM 6. Let G be a semi-simple Lie group and K a compact subgroup of G. Let X = G/K 
and X' = G' /K' , where G' is a subgroup of G and K' = G' fl K . Suppose a is the involution on G 
that fixes K and a' — o\x> ■ Then the natural embedding X' X is totally geodesic. 

The proof consists of noting that the hypotheses imply Cartan decompositions q = fi + p and 
q' = t' + p', and that p' is thus a Lie triple system in p: [p', [p',p']] C p, which is equivalent to the 
embdding being totally geodesic. 

Application of this theorem to the case at hand is immediate: Set G — SU (p, q), G' = SU (p' , q'), 
where p' < p and q' < q. Since any g' e G' satisfies g'*T p i^ q ig' = T p i^ q i it follows that the matrix 
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representation of g' has block-diagonal form, with a p' x p' block followed by a q' x q' . In order to 
obtain a corresponding element in G it is thus enough to insert a copy of between the 

two diagonal blocks of g' . In other words the embedding i : G —¥ G is simply 

( ^ \ 



.9 = 1 _ | ^9 

9q'xq' 



9p*xp' 

f-n—n'ri—n 1 



V g q'xq' J 

It is then clear that g satisfies g*T p , q g = T p _ q and is thus an element of G. K and K' , the 
maximal compact subgroups of G and G' , are afforded by the fixed points of the Cartan involution 
given by ac{X) = TXT, which agrees with cr* on su(p,q). That is to say, K = S(U(p) x U(q)), 
K> = S(U(p') X U(q')). 

With the above machinery in place, one is in a position to apply the results of Theorem 5 
directly, generating solutions for the harmonic map equation in this context. We carry out the 
procedure explicitly in the subsequent sections for N = 1 and a physically meaningful initial seed 
<7o- 

4.3 Einstein's vacuum equations: a 1-soliton calculation 

It is known [32, 21, 12] that the Einstein vacuum equations 

under the assumption of existence of two commuting non-null Killing fields K,K, reduce to the 
equations for an axially symmetric harmonic map / : A4 3 — > Af 2 . In the case one of the Killing 
fields is spacelike and the other timelike (outside a compact region in the spacetime,) the domain 
A4 is the Euclidean space M 3 . When both Killing fields are spacelike the domain is the Minkowski 
space M. = K 2 ' 1 , and in that case one usually calls / a wave map instead of a harmonic map (since 
the equations are hyperbolic and in fact can be written as a semilinear system of wave equations) . 
In either case the target Af is the real hyperbolic plane Hr, which is a symmetric space. By the 
results of the main theorem, the system is therefore completely intcgrable, and new solutions can 
be found from old ones by the vesture method [8] . This fact can be established directly, by taking 
G = SU(1, 1), K = U(l), and noticing H R = G/K. 

Einstein's vacuum equations in the stationary axisymmctric case can be written as a system 
for x and y, where for the spacetime metric g, x := g(K,K). Thus x > since K is assumed 
spacelike outside the axis of symmetry. Let b := *k * dK be the twist form for the Killing field 
K. It is not hard to see that the Ricci- flatness of g implies dh — 0, and thus we may let y 
denote a potential function for b = dy. Let (z, p, 9) denote Euclidean cylindrical coordinates on M. 3 . 
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Then the equations satisfied by x — x{z, p), y = y(z, p) are exactly those for an axially symmetric 
harmonic map into the upper half plane model of Hr. Indeed, the image of the map under the 
Cartan embedding of Hk into SX(2,R) is the following matrix 

1 I x 2 +y 2 y \ 

Now the Minkowski space R 3 ' 1 is clearly a stationary, axisymmetric solution of the vacuum 
Einstein equations and as such, it is a natural choice for a seed solution in the ISM scheme. It 
corresponds to the constant map 

x = l, y = 0, 

where we have let K denote the timclike Killing held Jj. The corresponding element in the sym- 
metric space G/K is therefore the identity matrix, 

go = «o(x) = / 2 x2 : M+ — ► G/K = SU(1, 1)/S(U(1) x U(l)). 

Note that any other asymptotically flat solution of the Einstein vacuum equations must have the 
same behavior as the above map, in the limit |x| — > oo. 

Let {TUk}^ =1 be a subset of C\K, and let Afc(x) be such that ro(Afc(x), x) = wt = w(Afc+jv(x),x). 
The set {Afc(x)}|^ 1 is the set of new poles, for a solution * of the Lax system, which is to be found 
using the vesture method: We have, trivially, \&ofc( x ) = ^2x2- From here one hnds the matrices 
Skj and the numbers ajy. One can then solve the linear system (50) to find the Ufc's in terms of 
the constants Wk and v/t, thus ending up with a 4iV-(real)parameter new solution q(x) as in (54). 
We'll do this now for N = 1: 

Fix a complex number w = is, s G M + and let Ai, A2 be the two roots of 

p(x, A) := A 2 - 2(z - is)X - p 2 = 0. 
Visibly, O = 2X 2 and by inspection, a simple corresponding solution ^ of the Lax system (9) is 

-^2x2- 

To solve (48), note that Skj = T and make the choice v\ — (a, S) 1 , so that v 2 — Tvi — (a, —S) 1 , 
where a, S are arbitrary complex parameters. It is possible to keep v\ free, noting there will be 
relations between the parameters (e.g., symmetries, quadratic constraint). For the purposes of this 
paper, the above assumption allows for simpler calculations. The resulting linear system AU* = B* 
is written as 







u\ 




—a 5 






u* 2 




—a —5 



(64) 
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Upon solving for each of the vectors u\,u 2 , the harmonic map is given by 

q(x) = I 2x2 - T-uivl - Y-U2«2- (65) 
Ai A 2 

Anticipating the Kerr solution as the outcome, we change from Weyl coordinates (p, z) to Boyer- 
Lindquist (oblate spheroidal) coordinates (r, 9) by setting 



p = \JJr — m) 2 + s 2 sin#, z = (r — m)cos9, (66) 

for real parameters m, s. Then the roots of p(x, A) are given by Ai >2 = [(r — m) ± is] [cos =p 1]. 
Solving the linear system (64) we obtain 

I ( a ( , ^) a( , n + B , ■ ) 

2D 1 5 c ^4 s n 5 / ^4 B \ 

\ ^sfcos^+l) r—m—is-' ^s(cos6— 1) r—m+is' 

where we have denoted 

B 2 A 2 



A:= |a| 2 -|<*| 2 > B := |a| 2 + |(5| 2 , D := 



4((r-m) 2 + s 2 ) 4s 2 sin 2 6' 
Plugging this into (65) we arrive at 

x + £ [-Asa8{iA{r - m) - Bs cos 0)] 



Q 



£ [450(5(1 A(r - m) + cos 0)] 1 + ^Elf! 

where F := A 2 ((r - m) 2 + s 2 ) - B 2 s 2 sin 2 9. 

The above (after adjusting its determinant to be one if necessary) is a harmonic map into the 
symmetric space Gi,i, for any value of the real parameters m, s and complex parameters a, S. Using 
the totally geodesic embedding of Qi i into SU(1, 1) we can view q(x) as an element of SU(1, 1) 
and consequently of SX(2,R), after a Cayley transform: Let 

1 i 
i 1 

Then q' := QqQ* e SX(2,R) for g e SJ7(1, 1), and thus, by (63), the Ernst potential £ = x + iy 
corresponding to a new stationary axisymmetric solution of vacuum equaions can be obtained by 
setting 

_ J_ A 2 ((r-m) 2 + s 2 )- B 2 s 2 sin 2 9 

X ~ ~q\ 2 ~ A 2 ((r - m) 2 + s 2 ) - B 2 s 2 sin 2 9 + 2s 2 (B 2 - A 2 ) + AsmA(r - m) - As 2 Bn 2 cos 9 

and 

q' 12 4sAn 2 (r — m) + ABnis 2 cos9 

V ~q 2 ~ 2 ~~ A 2 ((r -m) 2 + s 2 ) - B 2 s 2 sin 2 9 + 2s 2 (B 2 - A 2 ) + 4 sni A(r - m) - 4s 2 Bn 2 cos 9 

where we have set aS = n\ + in 2 . It is now evident that if we make the following identifications: 




B = a := \J m 2 + s 2 , n\ = — , n 2 = 0, 
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we would obtain 

r 2 — 2mr + a 2 cos 2 9 2ma cos 9 

r 2 + a 2 cos 2 9 ' ^ r 2 + a 2 cos 2 9 

which are precisely the expressions for the Kerr metric in Boyer-Lindquist coordinates. We have 

thus shown that give any two positive values m and s it is possible to choose the parameters a and 

5 in such a way as to recover the Kerr metric with total mass m and total angular momentum per 

unit mass a = Vm 2 + s 2 , namely it is enough to set 




Note that the Kerr solution thus obtained is necessarily naked, since a > m. 

It is likewise possible to obtain the Kerr-Newman solution to the Einstein-Maxwell equations as 
a 1-solitonic harmonic map into the complex hyperbolic plane He by dressing the trivial solution 
corresponding to the Minkowski metric, in the same manner as in the above, as will be demonstrated 
in the next section. 

4.4 The Einstein-Maxwell system 

It is quite remarkable that the Einstein-Maxwell system, 

R^-^g^ = K^; T^-^F^F^-g^F^F^- dF = 0; d*F = 0, 

which are the equations governing the interaction of the spacetime metric g with an electromagnetic 
field F permeating that spacetime, under the assumption of existence of two commuting Killing 
fields that also leave the field invariant, still reduce to the equations for an axially symmetric 
harmonic map / : A4 — > A/"', where M. is as in the above, and M' is the complex hyperbolic plane 
He [22, 12, 33]. Moreover, this target is a symmetric space, and indeed a non-compact Grassmann 
manifold £? 2 ,i = H c = G/K with G = SU(2,1), and K = S(U(2) x U(l)), so that once again, 
integrability is established and vesture method can be used to generate new solutions [6, 20]. 

To follow in the steps of the last example, it is computationally advantageous to make a change 
of basis and consider a unitarily equivalent representation of SU(2, 1). In particular, we use 

-i 

10- 

1 

The involutions are correspondingly defined as in (60), (61), so that the Lie algebra of G is given 

by 

S = su(2, l) = {Xe fll(3, C) | X*f + FX = 0}. 



G = SU(2, 1) = {g e GL(3, C) | g*Tg = f }, f := 
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We find a basis for g and confirm that it is eight-dimensional: Indeed, su(2, 1) = spanjJT 1 , . . . X 8 }, 
where 

X 1 = E 13 X 2 = E 31 X 3 = En — E 33 X 4 = i(En — 2E 2 2 + E 33 ) 

X 5 = E\2 + ii?23 X e = iE\2 + E23 X 7 = E 2 i — iE 32 X 8 = iE 2 i — E 32 - 

Here E^ denote the members of the standard basis for C 3x3 and C^- denote the structure constants 
of su(2,l), given in this basis by the commutation relations [X % ,Xi] = C^X k . The relations 
defining these constants are given in the table below. 





X 2 


X 3 


X 4 


X 5 


X 6 


X 7 


X s 


X 1 


X 3 


-2X 1 











-X 6 


X 5 


X 2 




2X 2 





-X s 


X 7 








X 3 









X 5 


X 6 


-X 7 


-X 8 


X 4 








3X 6 


-3X 5 


-3X 8 


3X 7 


X 5 




(anti-sym.) 






2X 1 


X 3 


X 4 


X 6 












X 4 


-X 3 


X 7 














2X 2 



Table 1: Commutation table for su(2, 1) structure constants 

We use Table 1 to easily determine the ±1 eigenspaces of a in q to be 

« = {Xeq \a«X = X} = sp'An{X 1 -X 2 ,X 4 ,X 5 -X 7 ,X 6 +X 8 } 
p = {Xeq \a„X = -X} = spsai{X 1 +X 2 ,X 3 ,X 5 +X 7 ,X 6 -X 8 }. 

It then follows (e.g. [7, p. 39]) that the Lie algebra q has the Cartan decomposition q — { © p, as 
described in section 2.5. 

To write down the Iwasawa Decomposition of q in our basis, we choose 

a = span {X 3 }, 

which is clearly a maximal subspace of p that is an abelian subalgebra of Q. For this a, we compute 
the root system to be as follows, where we have identified o* with R 1 , and ( ) denotes the linear 
span 

9 - 2 = (X 2 ), q-^ = (X 7 ,X 8 ), Q l = (X\X*), Q 2 a = (X^). 

Accordingly, A~ = {—2, —1, —1} and A+ = {2, 1, 1} (see section 2.6). The Iwasawa decomposition 
of the Lie algebra is then chosen to be q — n + ® a ® which lifts via the exponential map to a 
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decomposition for the group G = NAK, where K is the set of fixed points of a in G and A, N con- 
sisting of diagonal, unipotent elements, respectively, are the Lie groups obtained by exponentiating 
the algebras a, n + , respectively. 

We derive parameterizations for the subgroups N and A by exponentiating the corresponding 
Lie algebras. For the particular 4-dimensional representation of q that we have chosen, we have 

a = span{X 3 }, n+ = span {X 1 , X 5 , X e }. 

Therefore, an clement a e A can be parameterized, using one real parameter /i as 

aM=e"* 3 =diag(e",l,e-"), 

and similarly, n e N can be parameterized by one real parameter S and one complex parameter 
77 + 16: 



( r/ + i6 5 \ ( 1 r} + iO 5+Ur] 2 + 9 2 ) \ 



n(6,r,,0) = e sxl +* x '+ ex ' = cxp 



S 

it] + 6 




2 1 
iri + 6 

1 



/ 



We can thus find an explicit parametrization for the image of N A under the Cartan embedding, 
it will be of the form q — q(/i,S,r],9) = na 2 n* . This image is isomorphic to the symmetric space 
G/K, as claimed. For our purposes, however, it is useful to express this product in terms of the 
Ernst potentials 

v/2$ = f? + i0 £ = p> +iS = x + |$| 2 +iy. 



In particular, an element of the symmetric space is calculated to be 

' x + 2|<i>| 2 + \{y 2 + |<i>| 4 ) v^*(l ±\$\ 2 ) 

v^l + ^ + ^l 2 ) l + f|l>| 2 



P := ga(g)- 1 = 



X 



Notice that when $ = 0, one recovers an embedding of Q\,\ in ^2.1, as expected. We remark that 
the above representation of C/2,1 is unitarily equivalent to the one appearing in [13], p. 571, under 
the identification of our variables x, y, £ with <J/, $, H, Z, respectively. 

We would like to return this matrix to the original representation of SU{2, 1) in terms of V (as 
opposed to f), using conjugation by Q*. We thus obtain 



P := QPQ* = 



2|$| 2 -(£ + £) 



2\<S>\ 2 -\£\ 2 -1 2*(l-£) -i(\£\ 2 - £ + £ - 1) 

2$(l-£) -(2\<S>\ 2 +£ + £) -2*$(£ + l) 

i(|£| 2 +£-£-l) 2z$(£ + l) -(2|$| 2 + |£| 2 + 1) 



(67) 
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The role of P is to suggest which constants may be chosen when applying the inverse-scattering 
mechanism to generate solutions. 

We will now solve (48) for TV = 1 in this case, choosing the Minkowski seed qo(x) — hx3, and 
noting that Skj = T = diag(l, 1, — 1). Setting v\ = {a, ft, 7) , we have v<i = Tvi = (a, ft, — 7) , 
where a, ft, 7 are arbitrary complex parameters. The resulting linear system AU* — B* becomes 



1 



l«| 2 + 

m 2 + 



f + \S\ 2 ) -L- 



H 2 + 

M 2 + 



\ 2 + \s\ 2 ) 
\ 2 -\s\ 2 ) 



—a —ft 7 
—a —ft —7 



(68) 



A2 — Ai A2 — A2 

Our goal is to determine the appropriate choices for the complex constants using the matrix P 
derived above. We transform to Boyer-Lindquist coordinates to simplify the calculations, as before; 
note that the roots of Ai ; 2 of p(x, A) are exactly as in the previous example. It is easily verified 
that the linear system (68) has solution. 



/ 



a (s(cos8+l) 
R( iA 

A 



\ ^(s(cos6»+l) r 



r—m—is J 

B \ 
1 — m—is > 

B 



) P{ 
) 7( 



-iA 
sfcos 0—1) 

-iA 
s(cos — 1) 

-iA 
s(cos 0—1) 



+ 



r—m-\-is - 
B 

r—m-\-is ' 
B 



+is> ) 



where we have denoted 
A:= \a\ 2 + \ft\ 2 - 



B :-- 



|/3| 2 + l7| 2 



D := 



B 2 



A 2 



4((r-m) 2 + s 2 ) 4s 2 sin 2 B' 
Isolating each of the vectors Ui, 112, the harmonic map is given by 



<7(x) 



1 1 

^3x3 - —UlVl - T-U2V2 



Ai A 2 

j 1 + 8|a| 2 |7| 2 S 2 ^ [8a( g| 7 |2 s 2] 1 [_4 sa j( iA ( r _ m ) _ B S COS 6)] \ 

-p [-Asft~f{iA(r - m) - Bs cos 9)} 



V 



9(2,1) 
9(3,1) 



1 + 



8|l| 2 H 2 s 2 
F 

9(3,2) 



1 + 



/ 



Here, F := A 2 ((r - mf + s 2 ) - B 2 s 2 sin 2 6. 

We would like to show that the above six-parameter family of harmonic maps into C/ 2 ,i contains as 
a special case, the three-parameter family of Kerr- Newman metrics in Boyer-Lindquist coordinates 
(r, 9). To that end, we make the following identifications 



B 



-a := V 'm 2 + s 2 — e 2 . 



This gives rise to a system of equations for the constants a, ft and 7. Setting — n\ + in 2 and 
ftj — n 3 + in 4} one matches the 1-soliton solution generated from the Minkowski seed with the 
Kerr-Newman solution by choosing 



m 

n\ = — n 2 



n 3 



e 

2 " 4 



0. 
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One then obtains exactly the Ernst potentials for Kerr-Newman (see Eq. (21.26), p. 326 in [42]), 

* = : a £ = l : 3> 

r — la cos r — la cos 

for real parameters e,a,m. 

4.5 Higher-dimensional vacuum gravity and beyond 

The Qp^ q nonlinear sigma model described in the previous section has found applications in other 
settings (e.g., [40, 10]), and in particular in the study of higher-dimensional gravity. Moreover, 
explicating black- hole solutions in d-dimensional vacuum gravity for d > 4 has been of recent interest 
both in (minimal) supergravity and in string theory [19, 23]. In this context, stationary solutions 
possessing d — 3 rotational Killing fields have been extensively studied; imposing an additional 
timelike Killing field results in effectively two-dimensional theories, to which integrability techniques 
apply. 

For instance, the authors in [19] consider the Einstein vacuum equations for d = 5, having one 
timelike and two spacelike Killing fields, admitting a metric of the form 

ds 2 = g ab ( P , z)dx a dx b + e 2 ^ p ' z \dp 2 + dz 2 ), 

where detg = — p 2 . The Einstein equations then reduce to 

d p U + d z V = 

U := (pdpgg- 1 ) , V := (fagg- 1 ) 

d p u = -^ + ^tv(U 2 -V 2 ) 1 d lV =^-tr(UV). 
2p Hp 4p 

The first two equations of the system comprise a principal chiral field model into GL(3, K), coupled 
to equations in v(p, z) which may be solved by quadrature once g is determined. The Zaharov- 
Bclinski technique is employed to produce Kerr, Myers-Perry and black ring solutions. On the 
other hand, one could look for a sigma-model representation of the field equations (that is to say, 
a harmonic map into the symmetric space SL(3,M.)/SO{2, 1)) and employ the solution-generating 
methods of this paper to address the case of a non-diagonal initial seed. This and other applications 
will be pursued in a forthcoming paper. 

Further progress on this topic may also result from extending the integrability results of this pa- 
per to situations where the target is not a symmetric space but a homogeneous space, in supergravity 
descriptions of M-theory, for example [15, 16]. These ideas will be pursued elsewhere. 
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4.6 Summary & Outlook 

By establishing the integrability of harmonic maps from effectively 2-dimensional domains into 
Ricmannian symmetric spaces G/K for real semisimple Lie groups G, and showing how the dressing 
technique can be used to generate new solutions from known ones, our paper goes beyond the 
current literature in providing a general framework for the study of harmonic maps into noncompact 
symmetric spaces commonly appearing in mathematical physics. Examples we consider suggest that 
this approach may also make generating solutions for other effectively two-dimension geometric 
field theories more tractable. Future directions for research on this topic include the possibility of 
dressing with poles on the real line, which would give rise to a different class of singularities for the 
dressed solution (e.g. black holes vs. naked singularities), and extending the integrability results 
of this paper to more general targets, e.g. homogeneous spaces. 
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